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^ Abstract 

We approximately compute the bispectrum induced on the CMB temperature by fluctuations in the 
CN standard recombination epoch. Of ah the second order sources that can induce non-Gaussianity 
^ during recombination, we concentrate on those proportional to the perturbation in the free electron 

density, which is about a factor of 5 larger than the other first order perturbations. This term 
I induces some non-Gaussianity by delaying the time of recombination and by changing the photon 
2 diffusion scale. We find that the signal is not scale invariant, peaked on squeezed triangles with 
CZ3 the smaller multipole around the scale of the first acoustic peak, and that its size corresponds to 
I— —I an effective /nl — —3.5, which could be marginally detected by Planck if both temperature and 
polarization are measured. 
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\o 1 Introduction 
cn 

CN Non-Gaussianities in the Cosmic Microwave Background (CMB) have become very important in the 
last few years. On the observational side, there has been a huge improvement in the sensitivity with 

<^ experiments such as WMAP. Upcoming experiments such as Planck will be sensitive to even smaller 
t> non-Gaussian signals. On the theoretical side, standard slow-roll inflation predicts an extremely 

^ small amount of non-Gaussianities [1], undetectable by next generation experiments. However, 
^ several large modifications to the standard inflationary picture have been proposed both in the 
case of single field infiation jH El SI |5] , where all the models have been unified in an effective field 
theory description in [6! , and in the case of multi-field infiation [TJ |H] . Recently even a consistent 
bouncing cosmology has been proposed [9l [TOl [H], which, though significantly less compelling than 
inflation, does predict a possibly detectable level of non-Gaussianities [Tl] . A detection of primordial 
non-Gaussianities would therefore imply a radical departure from the standard and familiar slow- 
roll inflation scenario. Non-Gaussianities are probes of the interactions of the inflaton, and therefore 
contain an unprecedented amount of information on inflation. If detected, we would have to abandon 
the standard slow-roll inflation picture, but would be left with new information to understand the 
dynamics of the inflaton. 

Of course, non-Gaussianities in the CMB are not generated only by the inflaton, there are 
many contributions coming from the post-inflationary evolution. These can be roughly subdivided 
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into effects wliicli are important around or before recombination, and foreground contamination or 
secondary anisotropies generated after recombination. The latter include scattering secondaries, 
such as the thermal Sunyaev-Zel'dovich (SZ) effect [KJ and the kinetic SZ or Ostriker-Vishniac 
(OV) effect [U]; and gravitational secondaries, such as weak lensing and the Rees-Sciama (R-S) 
effect |14] . In order to be able to fully exploit upcoming CMB experiments, it is necessary to 
understand quantitatively the contributions from all these post-inflationary mechanisms. 

The effect of the secondaries on the CMB bispectrum has been already extensively analyzed in 
the literature. The thermal SZ effect and its frequency dependence are investigated in [15]. The 
R-S effect is studied in [16 . Gravitational lensing and the ISW effect are considered in [T7] and 
|18] . The OV effect is analyzed in e.g. p9] and pO]. The signal-to-noise for upcoming experiment 
generated by these secondaries is very large and it will have to be taken into account when looking 
for primordial non-Gaussianities. 

Let us now consider non-Gaussianities generated from around recombination and concentrate 
on the 3-point function or bispectrum. Due to translational invariance, the bispectrum in Fourier 
space is a function of three wavevectors which must add up to zero forming a closed triangle. For the 
squeezed triangle limit, where one of the modes is well outside the horizon, the bispectrum generated 
around recombination has been obtained exactly in |21j|. Oversimplifying, we can imagine that the 
non-Gaussianities are given by the following relation between the observed curvature perturbation 
C and a gaussian random variable (g-. 

C(x) = Ux) - If^l- iUx)' - (C,(x)^)) , (1) 

where is a measure of the level of non-Gaussianity, and the superscript loc, which stands 
for local, refers to the above local-in-space structure of the non-Gaussiantiy. Then, in [21^ it was 
found that the non-Gaussianity in these squeezed triangles corresponds to an f^^' < 1, which is 
undetectable even by Planck. 

The purpose of this paper is to go beyond this simplified limit, and take a step further in the 
computation of the non-Gaussianities generated around recombination for all triangles. The full 
calculation of the bispecturm is a very hard task. In order to obtain it, one is forced to solve the 
second order Boltzmann and Einstein equations, which are very complicated. We give those in a 
companion paper [22], following and correcting the equations given in [23] and [2l], and including 
the effects of the free electron number density perturbation. This we calculated at first order 
in the approximation of the Peebles' effective 3-level atom [25]. In |22] we showed that around 
recombination, the amplitude of the perturbations to the free electron density 5e is enhanced by a 
factor of ~ 5 relative to the baryon density perturbations due to the relatively small timescale of 
recombination Q 

Without accounting for the enhancement of 5e, one very naively expects that the second order 
evolution will lead to a CMB bispectrum with an ~ 1, which for an experiment like Planck 
corresponds to a signal-to- noise ratio of ~ 0.3 (see for example [26]). Consequently, we may expect 
that the enhanced 6e will lead to an observable three-point signature in the CMB corresponding to 
an f^^' ~ 5. In this paper, we therefore concentrate only on those non-linear perturbations which 
are induced by 6e- 

^As usual, we denote the change in a given quantity q by 5q, while the fractional change we denote by 
5q = Sq/q 
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Notice that, as long as the tight-coupHng between photons and baryons holds, the CMB tem- 
perature does not depend on the electron number density Tie, and the contribution from 6e to the 
second order CMB anisotropics vanishes. This means that the effect of 5e is most important during 
recombination. 

At first order, the solution for the CMB temperature anisotropics O = 6T/T, once expanded in 
multipoles, is given as an integral over the first order source S^^^ and the photon visibility function 
g{f]) |Z3- I^oi' ^ mode of wave vector k, the solution reads: 

Gp)(A;,7?o)~ / drjg{rj)S^'\k,r^)ji[k{r^o-ri)] , (2) 
Jo 

where ji is the spherical Bessel function, r] is conformal time and ryo corresponds to r] of today. The 
order in perturbation theory is given as a superscript in parenthesis and will be dropped whenever 
that does not cause confusion. If S^^^ is expressed using the photon monopole, dipole, and quadrupole 
without accounting for photon diffusion, then including the effects of photon diffusion approximately 
amounts to replacing S^^^ — > S^^^ exp[—k'^ /k'jj], where k£) is the photon diffusion scale (see for 
example Section 8.5 of |28j). 

The electron density rie enters only in the visibility function and the diffusion scale, which means 
that 5e multiplies the first order source. Thus, any perturbations to rig will affect Q at second order 
only. This can be understood also intuitively. In a homogeneous universe, before, during and after 
recombination the radiation temperature decreases as a~^, o being the scale factor, irrespective of 
the electron density. This means that a perturbation to the electron density changes the position 
of the last scaterring surface and the mean free path of the photons before recombination, but not 
the observed radiation temperature. At first order we perturb Ue and keep the other quantities 
unperturbed, and therefore the CMB anisotropics are not affected by 6e at this order. 

From ([2]), we expect that the expression for the second order temperature anisotropics after 
including the perturbation to rif, will be schematically given by[^ 



eS(fc,r?o)~2 rdw(r?) 
Jo 



k^ 



s^'\k,7^)mm-v)] , (3) 



where 5g and are the fractional perturbations to the visibility function and diffusion scale. When 
later in the paper we perform the full calculation, we will see that the above schematic equation is 
not far from the right expression. 

The paper is organized as follows. In sec. [2] we elaborate on eq. (|3| and provide simple estimates 
of the induced bispectrum. Then, in sections [3] and |4] we go on and present the full calculation of 
the bispectrum concentrating on modes inside the horizon so that we expect to be able to neglect 
second order metric perturbations. We summarize the result of the analytical calculation in sec. [5} 
In sec. [6] we show the results of the numerical integration and describe the induced signal. We 
conclude in sec. [7j In App. |A]we summarize some useful formulas about rotation invariance. 



^In the paper we expand the perturbations as 



/ (x , p\ 77) - /(°) (p\ 77) + /( 1) (f , p\ ,7) + i / (2) (f , ^) 
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2 (Not so) simple estimates 
2.1 Mechanisms generating 

Let us begin to analyze the mechanisms by which perturbations in the number of free electrons can 
induce a bispectrum in the CMB. We restrict our analysis to modes which are faster than the horizon 
scale at recombination (k > 0.005 Mpc~^ ), and slower than the scale associated to the width of 
the visibility function (k < 0.27 Mpc~^ ). This is the region from where we expect that a large 
three point function can be generated. For modes which are out of the horizon at recombination, 
the calculation of [21] which is valid in the limit in which one of the modes is outside of the horizon 
and the other two modes are much faster than the first one, shows that the expected signal is 
equivalent to < 1, which is very small. Similarly, we do not expect a large signal coming from 
configurations where all of the modes are out of the horizon and of comparable size. We do not 
expect a high signal also from modes with k > 0.27 Mpc~^ because in this case the mode is faster 
than the timescale of recombination, and there is effectively an averaging out of the perturbation 
|22] . Furthermore, these very high k modes lie outside the reach of experiments such as Planck. 

As eq. ^ suggests, 6e will affect the CMB bispectrum in three different ways: by perturbing 
the physical time at which recombination takes place, by perturbing the photon diffusion scale, 
and by perturbing the probability for the CMB photons to originate from different times within 
the recombination era. This last mechanism could be relevant because diffusion damping during 
recombination smoothes out the perturbation associated with a given mode on a time scale corre- 
sponding to approximately the period of oscillation of the mode itself. The later in time photons 
originate from (even if an order one Hubble time after the peak of the visibility function) , the more 
the anisotropy is suppressed. We will find that this effect is subdominant. 

The different contributions can be extracted by writing eq.s ^ and ^ as 

rvo 
Jo 

xS^^y{ri + 6r]r{ke))ji[ks{m - V - '^??r(4))] , (4) 

where 5"^^^' is the first order source without the exponential damping; ne is the free electron number 
density; Cg is the sound speed; A;^ is the wavenumber of the source; kf, is the wavenumber of 5e'-, g^^^ 
is the unperturbed visibility function; and a convolution such that kg + kg = k is implicit. 5r]r{ke) 
is the perturbation to the position of the peak r]r of the full visibility function 

g^''Hv)+Sg{v,ke) , (5) 

where 6g is the perturbation to the visibility function. Still in eq. Q, we have split the perturbation 
to g into two pieces: 

Sg = -g{'n)^'nr{ke) + 6g{ks, 4) • (6) 

The first piece gives approximately the time shift of the (otherwise unperturbed) visibility function, 
while 5g takes into account the modification of its shape. 

By shifting the variable of integration, in eq. ^ we have moved Srj^ inside the first order source. 
In the limit in which the visibility function can be roughly thought as a 5- function centered at r]r, 
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the effect of 5rjr is to evaluate ji and all the time dependent quantities in the first order source term 
at a shifted time. After this time shift has been taken into account, another possibly important 
effect of the perturbation of g is its change in shape. As we will later see, this effect will turn out 
to be unimportant, which we will show in the estimates that follow. Even in the case in which a 6e 
mode does not change the position of the peak of the visibility function at recombination, it can still 
distort its shape, by for example making the width smaller. Because in the time scale of the order of 
a few widths of the visibility function the source decays relevantly (even for relatively low ks), the 
actual location around recombination from where most of the CMB photons originate, is relevant in 
determining the magnitude of the CMB anisotropy. If there is more probability in the central part 
of recombination, then the anisotropy is larger, while in the opposite case, it is smaller (see Fig. [T]). 
Clearly, this effect vanishes in the limit where the diffusion damping scale goes to infinity, because 
in this case the size of the source is independent of the place from where the photon originates. It 
also vanishes in the limit where the 6e mode is oscillating on a time scale either much faster than 
the diffusion damping scale of the source mode, because in this case an averaging occurs, or much 
slower, because in this case it will just shift the peak of the visibility function, but not change its 
shape. Of course, for a generic 6e perturbation, a part of the above effect is taken into account by 
the shift in the position of the peak r]r of the visibility function. In order to take into account the 
effect due to the change in its shape, we define a perturbation to the 'effective Area' of the visibility 
function after the time shift has been subtracted: 



where, as explained, the weighting by the diffusion damping takes approximately into account the 
importance of the location of the probability to create a CMB photon, for the resulting size of the 
CMB anisotropy. 

Before we procede we would like to have an expression for the perturbation to the visibility 
function, 5g, wich includes the oscillatory character of 5e- This can be taken into account if one 
perturbs 6g not simply with the amplitude, 5e, but with a multipole expansion of 6e{ke,r]) exp{ike-x). 
The angle-independent part of the multipole expansion of exp(i/ce • n(r/o — v)) is oc jz(fee(??o — ^0)- 
We therefore use the following effective perturbation to g: 



ji,ff{les)Sg{ke,'n) = g^°\r]) [6e{ke,'n)ji^g{ke{rjo - rj)) - dr]'f{r]')6e{ke,'n')Xg{ke{r]o - v')) ] , (8) 



with /eff = ^e(% — ??r)) where we have pertured g with Se{ke,ri)ji^g(ke{r]o — rj)). This choice of 
perturbation approximates the result of our full treatment (cf. eq. (|59])). 

The total perturbation to g at late times is large because there is still quite a bit of optical depth 
after recombination (r = 0.04 between rj = 360Mpc^^ and r]o) and de grows once it starts falling in 
the CDM potential wells [22]. This could naively give a large effect. However, in order to accumulate 
a good fraction of this optical depth the photons need to travel a large distance comparable to rjQ. 
This is possible only for k being to a very good approximation perpendicular n, since otherwise the 
oscillating 6e averages out along the line of sight. This implies that for a fixed \k\ and I, there are 
very few k directions which contribute to the signal. This effect is taken into account in eq. ^ by 
the spherical Bessel function. This suppression of 6g at late times renders the (5Area effect small as 
we will see below. 



(5 Area 



j drig^^\vi)5g{r]) exp[-A;^/fcg)(r?)] 
/d7?5(0)(r/)exp[-fc2/A:2^(r?)] 



(7) 
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Figure 1: Schematic representation of the effect of a perturbation 5g to the shape of the visibihty 
function for the case high-Zcg and low-A;^. Notice that the location of the peak is unchanged, so that 
Sijr = 0. However, g is much more peaked around the central value than in the homogenous case. 
This means that more CMB photons originate from the region near the peak, where the power of 
the mode is not yet damped (in red we plot the diffusion damping). This means that the resulting 
CMB anisotropy is larger. The shaded black region of g^^^ represent how much the 'Effective Area' 
in the homogeneous case differs from one and it is a measure of how many CMB photons originate 
from the region where the mode is already beginning to be damped. We see that g^^^ + 6g has a 
larger 'Effective Area', implying a positive (5 Area and the possibility for a larger CMB anisotropy. 



From eq. Q, we obtain an approximate expression analogous to eq. ^ that we can investigate 



semianalytically and that explicitly differentiates between the different mechanisms generating @ 



(2). 



Vo 



d'i]g{T]) exp[-k'^Jkl,{rj)] 



D + '^jt^h + (5Area + 5rir <9^ 



(9) 



where 5kg = 6r]rkD/k£) is the perturbation to ko due to the timeshift Jr/^ in the visibility function 
which makes ku to be evaluated at a slightly different time than in the unperturbed universe fnotice 
that k£){vi) oc sjueirf) ), while measures the perturbation to /cd due to 5ne directly]^ It is 
quite straightforward to realize that, in the limit in which the timescale over which the source varies 
is much longer than the width of the visibility function, the above equation gives an approximate 
expression for the dependence of on the JArea. This expression is expected to be correct up 

to order one corrections, which is enough for the purposes of the present section. As we will later 



show (see later in eq. (12|), the 3-point function will be sensitive to the variation of the angular 
power spectrum Ci with respect to the quantities perturbed by 5e- In the case of the phase shift (5t/j., 
this gives rise to a term proportional to {—Cska5r]r)s\n{csksr]r), which, once integrated in Fourier 
space to obtain the Cfs, gives rise to approximately an effect of order (—CsksSrjr) /3. In the regime 
of interest, the factor of ~ 1/3 takes approximately into account of the difference from integrating 
over ks the source expression that now contains cos{csksr]r) sm(csksrir) instead of cos'^{csksr]r)i in 



^Notice that we denote by Sku the full perturbation to k^, due to both of the effects, and not just to Srie- 
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both cases multiplying other kg dependent terms. We verified that the perturbation to r] entering 
in ji leads to a negligible contribution, and we therefore drop it here. 




Figure 2: The first-order perturbations responsible for the bispectrum generated from recombina- 
tion. The computation was done in synchronous gauge, and the normalization of the perturbations 
is the one in CMBFAST [29]- on super horizon scales C = 1- represents the perturbation to the 
diffusion scale at the peak of the visibility function (77 ~ 288 Mpc) due to drie] 5kg is the perturba- 
tion to ki) due to 6r]r - the shift in the position of the last scattering surface. (5 Area approximately 
represents the perturbation to the probability that a photon originates from the last scattering sur- 
face before the perturbation decays due to diffusion damping. The timeshift Srjr also gives rise to 
a change in the phase at which we evaluate the first order source, which is schematically given by 
{—Csks6r]r/3). The plot is for kgijQ = 3000. From the plot we can see that the largest contribution 
to the second order temperature anisotropics for the given ks is from kgrjo ~ 200, i.e. from around 
the first acoustic peak. Note that the two largest contributions to - the perturbations oi k£) 
coming from 5e and the change in the position of the last scattering surface, partially cancel each 
other. 

Using our results from |22] for 5e, in Figures |2] and [s] we plot the above four first-order quantities. 
All calculations were done in the synchronous gauge. 

When we compute the bispectrum, we compute the correlation function of three modes. In 
order for this to be non zero, one of the modes has to be taken at second order, obtaining an 
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Figure 3: The same as in Fig. [2] but for k^rjo = 200. As one can see, all effects for high suffer 
from some kind of suppression as described in the text. 

expression of the form (neglecting multipole indices for simplicity) 

(e(2)(fci)e«(fc2)e(i)(fc3)) ~ (<^«(4)e«(fc,)e«(fc2)e«(fc3)> , (lo) 

where ke + kg = ki. As shown in Fig. |4j when computing the expectation value, we pair each first 
order perturbation with the others. This forces k2 = —kg and k^ = —ks, while the sum of the three 
momenta ki, /c2 and k^ must form a closed triangle by translation invariance. As we can see from 
Fig. [2] for a given high ks, the signal grows as we make ke lower |^ In the opposite case, for a low 
ks, we see in Fig. [3] that the signal never grows as we make kg larger. This makes us suspect that 
most of the signal comes from combining a low kg with a first order source term at high ks- This is 
actually true, as we will confirm in the full calculation: the signal is peaked on squeezed triangles 
as the one shown in Fig. |4j 

The two largest contributions to the combination of a low ke and a high kg mode (see Fig. [2]) are 
given by the perturbations to the diffusion scale, 5^^ and 5^^, which are comparable in magnitude 
but out of phase. These quantities are not enhanced by themselves, but they are multiplied by a 
factor of 2k1/k'j~f in the expression for the second order temperature anisotropies Q, which boosts 

"^This is true only until he is within the horizon. For fcg outside of the horizon 5e becomes irrelevant. In 
this paper we refer to low-fc modes as modes that project on the sky as multipoles with I ~ 200, while by 
high- A: modes we mean modes that give rise to multipoles of order thousands. Similarly by a low-/ we mean 
an I ~ 200 and by high-Z an I ^ thousands. 



8 




Figure 4: When computing the bispectrum, we take the expectation value of three G modes, and 
one of them has to be taken at second order in order not to have a null result. Each of the two 
first order perturbations contained in the second order mode, approximately 5e and Q^^\ need to be 
matched with one of the two first order G modes, forcing ke = —k2 and kg = —k^, where ki = k^ + kg. 
The sum of /ci, k2 and k^ must be equal to zero, so that the three wave vectors form a closed triangle. 
The same is done for all the symmetric combinations. 



their individual effects on the bispectrum, although their net effect is nearly completely cancelled. 
This cancellation can be easily understood since Ue is larger in overdensities. This means that in 
the presence of a \ow-ke overdensity, the high k modes (implying high ks) will be less supressed 
due to the reduced diffusion scale (i.e. positive (5^^), but they will have more time to decay (i.e. 
negative 5kg) until recombination takes place. This can be understood also in the following way. In 
the limit in which the 5e mode is much slower than the timescale of recombination {k < 0.1 Mpc~^), 
we expect that recombination happens in the same way as in the unperturbed universe, just a bit 
time translated. For this reason, around recombination, rie will be just the time translation of its 
unperturbed value. Since k'^^ ~ aTne{'r])a{rj)/r], where or is the Thomson cross section, after having 
taken into account of the time translation of ng, the piece proportional to 6ne disappears, and one is 
left only with a perturbation to the scale factor and to tj, which give a very small effect. Notice that 
in this same regime, Ue is determined by the local value of (or equivalently of the temperature) 
with the same relationship as in an unperturbed universe: 6n(, ~ dnbhe/fib- As one can see in 
Fig. 5 of [22], this is in fact the case even for modes very well inside the horizon, implying that the 
cancellation between and begins to be milder only for relatively high I modes. 

This is not however the whole effect in the regime of low ke and high kg. Still from Fig. [2| we 
expect that for this combination of /c's, the phase shift due to 6rjr in the source should also play an 
important role for the bispectrum. As we will later show, the full calculation involves an integral 
in time, which suppresses the effect from the phase shift with respect to the naive expectation from 
Fig. [2] due to the fact that in this case the integrand changes sign. Finally, in this regime, the 
perturbation to the area is not significant, which tells us that, as expected, for low kg, 5g is well 
approximated by just a timeshift in o. 

(2) r~i 

The contributions to G^ in eq. ( 9 ) from the combination of a low kg source mode and a high ke 
be are shown in Fig. p3l Since kg is sma 



I, the factor 2kg /kj-, supresses the contribution of (5^^ and 5kg 
in this case. Also the time shift is quite small. The effect that comes from the perturbation of the 
shape of the visibility function and therefore to the probability for the CMB photons to originate 
from different times within the recombination era is also shown in Fig. [3j This is what we called 
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6 Area and described at the beginning of the section. This effect does not become important even for 
very high ke that change appreciably within the width of the visibility function, due to cancellation 
effects as described around eq. ([S]). 

For a high kg perturbation of order 1, we can estimate the effect it will have on the bispectrum. 
Starting from eq. (10), we can write the ratio of the induced bispectra in the two different {ke,ks) 
limits as 

(perturbation(A;high) ~ l)©(Mow) X e(A:high) X 0(A;iow) 1 „^ 

p ~ r T ~ U.S , (llj 

2^4o(A;iow)e(fchigh) X e(/chigh) X e(A;iow) 13 Exp [-k^^Jkl 

where A;high% — 3000 and kio„r]o ~ 200, and where we took the value of S^^ from Figure [2j Notice 
that the reason why the above ratio is order one is not because the high-Zce perturbation is large, but 
because the ratio involves 6(A;iow)/0(fchigh) which is enhanced by the exponential damping. This 
means that not-enhanced first order perturbations that are not exponentially suppressed at high k, 
might give rise to an enhanced bispectrum. All of our high-Zce perturbations are suppressed, as is 
the case for (5 Area due to cancellation effects, and therefore we do not expect a large signal to noise 
from this regime. 



2.2 Estimates for the bispectrum 



As discussed in Section 2.1 , we expect a possible enhancement to the bispectrum from recombination 
only for modes with / between approximately 100 and 3000. A simple approximate formula for the 

l~l (2) 

3-point function can be obtained by noticing that expression ( 9 ) for the second order is simply 
given by the product of a perturbation to a quantity 'X' induced by 5e times the derivative of the 
first order source with respect to X. Therefore, in this case, in the limit in which Sx does not 
vary appreciably during the width of the visibility function, we can approximately write the 3-point 
function of the multipoles of the CMB anisotropy a^^'s as 



{ai 



3^3/ 



->mim2mz -h+h+h 



where 



1 dC, 



2dlnX 



+ 5 perm. , 



(12) 



(13) 
's (see 



^hhP^"^'^ is Gaunt integral, and the factor i'l+'z+'s comes from our choice of phase for the 
eq. (15) below). In the flat sky limit the Gaunt integral just reproduces (27r)^J^(Zi + I2 + h)- We 
will verify explicitly this formula when later we perform the full calculation. In the above equation 
we defined 



in analogy with 



Sx. 



air. 



Im 



2^2 



27r2 



d'kYiUk)ak)6xik,r]r)jiikr]o) 



d'kY;^{k)C{k)S^'\k,Vr)ji{kvo) 



(14) 



(15) 



where rjr denotes ij at recombination. Here we have introduce the primordial curvature perturbation 
C{k). In this case the first order quantities like 6{k, rjr) or S^^\k, rir) are meant as transfer functions. 
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We can perform the angular integration in the expectation value in (13) to obtain 



dkk^P{k)ei,{k,r]o)ji^{kr]o)5x{k,Vr) x 7, 



1 dCi, 



2 

vr 



dkk^P{k)Qi^{k,rjo)ji2ikr]o)6xik,r]r) x 



2ainX 
1 dCi, 
2dlnX 



+ 



(16) 



+ 2 perm. . 



When later we perform the full calculation, we will see that the above equation is a good approx- 
imation in the limit in which 6x does not vary within the width of the visibility function, which 
corresponds to multipoles of order a few thousands. 

The above expression can be evaluated numerically with the approximate expression for ji given 
in eq. 9.71 of [30j, and using the approximation for the monopole and dipole of the first order source, 
given in eq. 8.24 and 8.26 of [28] (where we drop the integral contributions for simplicity). We also 
include the exponential damping. The sign of the reduced bispectrum is chosen such that a positive 
bispectrum in the local model implies a negative /^l' ™ the Sachs- Wolfe regime, and vice versa. 
Of course this expression will not be very accurate, due to the approximate treatment of diffusion 
damping in the first order source, to the additional averaging across the visibility function, which 



will reduce the effect from Srj^ and 5 Area, (see eq. (59)), and to the approximation in the dependence 
of 6(2) on 5Area. Stih, we will see that it reproduces quite well the results of the main calculation. 

For the perturbations due to recombination we are dealing with, as we saw in eq. ([9|, 6x runs 
over {6ri entering in S'^^)'}, 5^^, 5kg and (5Area. Since later when we do the full calculation, we find 
that most of the signal is peaked on squeezed triangles with sides li ~ 200 and /2 ~ ^3 ~ 3000, we 
concentrate our discussion directly on this limit. The spherical Bessel function ji in the integral in k 



in each term of (16) forces the k of 6x{k) to be approximately equal to l/rjQ. This means that when 
a term has I = h, as the first one shown in (16), it represents an effect due to the low- A; part of the 
perturbation 6e, while when it has / = as the second one shown in (16), it represents an effect 
due to the high-A; part of the perturbation 6e- For each of these separate terms, in Figures 5][7 we 



plot the approximate contributions to the signal-to-noise per triangle of the three-point function in 
the squeezed limit as obtained from (16). The plots are for li = 200. 



Let us start discussing the contribution arising from a low kf, and high kg. The contributions 
from {dkf) due to Sn^} and 6g are shown in Fig. |5] As discussed before, the low-k^ contributions 
nearly cancel each other, which is reminiscent of the fact that for this low-fcg modes recombination 
happens as in the unperturbed universe, just a bit time translated. This can be also confirmed by 
comparing the bispectra they generate: the effect from the low-fce S^g is shown separately in Fig. 6 
and one can see that it is very nearly equal to minus the effect of the low-fee (5^^ shown in Fig. 5 
Since ke is fixed and small, we have kg ~ h/vo- Notice also that from eq. ^ we could have expectec 
that each of these two contributions should scale as I2 k^, which is nicely confirmed by these plots. 

From Fig. |6] we also confirm that the phase shift of the source due to 5r]r has a contribution to 
the bispectrum which increases linearly with ks and is non-negligible, as suggested by our analysis 

r~i 

of its contribution to (see Fig. 2). In the full calculation this linear growth will become milder 
at higher /c^'s because there is a time integral which tends to cancel out the effect. 

Let us now consider the other interesting regime. The combination of a fixed low kgijQ = 200 
with a high kg is shown in Fig. [Tj As anticipated (see Fig. |3]) we find that there are no important 
contributions in this regime. 

Let us now explain the signs of the contributions proportional to dk^. We can concentrate 
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Figure 5: Approximate contributions to the signal-to- noise of the three point function per triangle 



in the squeezed limit as obtained from (16) for li = 200 as a function of I2 = h- The contribution 



to the bispectrum from 5/^^ comes from the perturbation due to 6ne- The different contributions to 
the bispectrum from 6g are shown in Figures |6] and [7j 
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Figure 6: The different contributions (obtained from the approximation given in eq. (|16|)) to the 
signal-to-noise of the bispectrum per triangle due to a low- A: 6g for li = 200 as a function of h = h- 
The sum of all these curves is equal to the \ow-k 5g curve in Fig. |5j Due to the ensemble averaging 
of the primordial fluctuations, we have ke ~ h/voi ^^^d ~ h/vo- The 5ko due to 6r] is comparable 
and opposite in sign to the dko due to 6ne (see Fig. [5]). The scalings of the bispectra due to the two 
contributions to Ska follows kg oc ^ implied by eq. Q. The change in the phase of the first order 
source due to the shift of the position of the last scattering surface brings an oscillating contribution 
proportional to kg, whose amplitude grows linearly with kg, as expected from eq. ([9]). 
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Figure 7: The different contributions (obtained from the approximation given in eq. (|16|)) to the 
signal-to-noise of the bispectrum per triangle due to a high- A; 5e for li = 200 as a function of h = h- 
Here we have ke ~ h/vo^ and ks ~ h/vo- Although a high-Zce perturbation of order 1 can produce a 



large bispectrum (cf. eq. (11 )), all high k^ effects due to recombination are suppressed and therefore, 



we expect that the bispectrum generated in this regime to be small. 
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on a low-Zce perturbation since this is the regime in which they are relevant. In this hmit, the 
most important terms of the bispectrum in (16) are just the ones where k^, ^ h/rjo, which we can 
now schematicahy write as 0;^ x ^jOzjl • This tehs us how the amphtude of the high-frequency 
modulations changes in the presence of a long-wavelength mode. If the bispectrum is positive, high- 
l modes have an enhanced amplitude in hotter long-wavelength regions of the sky, while they have 
a smaller amplitude in colder long- wavelength regions P) 

(2) 

With this intuition in mind, let us now check whether the different mechanisms generating 
produce bispectra with the signs obtained with our approximate equation ( 16 ). Let us concentrate on 
the interesting regime, corresponding to the low-fcg (correspondoing to low li) 5kg terms. Notice that, 
since G;j(ryo) ~ S^^\k = ^i/r/o,?7*), a low-Z mode at the scale of the first acoustic peak corresponds 
to a positive temperature perturbation on the sky, which is also associated to overdensities. But as 
we discussed in our former paper ^2] , in overdensities rig is larger and thus is positive around the 
same scale (as can be seen by looking at the phase shift term in Fig. |2|. Therefore, in the presence 
of such a low-/ mode, the photons have more time to diffuse, which decreases the amplitude of the 
high-/2 Sza- This means that we can expect (5^^ entering in (i.e. the 5kg effect) to generate a 
negative bispectrum, confirming what we find in Fig. [6} This conclusion holds even for higher li 
(around the second acoustic peak), since there 5^^^ and Q^^^ have again the same signs. 

In the case of the 5kj^ perturbation generated by Jrie, the perturbation is clearly positive in low- A; 
overdensities. By exactly following the same argument we made for 5kg, it is immediate to see that 
the induced bispectrum in this case is positive, confirming what found in Fig. [5] 

From these estimates we see that the bispectrum produced by perturbations to the recombination 
history can be quite large, though some cancellations tend to reduce its net size. This clearly deserves 
a more accurate computation, which we are going to do in the remaining of this paper. In order to 
do this, we will work well inside the horizon, so that we expect to be able to neglect second order 
metric perturbations; and we will also drop, in the second order photon Boltzmann equation, all 
those quadratic second order terms that are not proportional to 5e- This procedure will grasp the 
leading effects proportional to 5e which are expected to be enhanced. 



3 Second order temperature anisotropies 
3.1 Line of sight integral 

Let us procede to find the exact expression corresponding to our intuitive expression ([s]). We follow 
the conventions of [22]. As explained, we work at first order in the metric perturbations, restricting 
our calculation to be valid for scales much shorter than the horizon at recombination. In this 
analytical part we work in Newtonian gauge. The line element is given by: 

ds"^ = a^{ri) [-(1 + 2^)ci7?2 + (1 - 2^)5ijdx'dx^] , (17) 

where we ignore primordial vector or tensor modes. 

^This means that a positive bispectrum implies a positive skewness of the one-point distribution function 
of the CMB temperature anisotropies. 
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We start with the second order brightness equation given by eq. (76) in ^22j. It is an equation 
for the second order perturbations A^^) of the photon energy density: 



(18) 



for i = 1,2 , where 



F{x,p\rj) = F^^^\x,p\T]) + F^^\x,p\r]) + ^F^^\x,p\r]) 



is the gauge invariant photon one-particle distribution function written to second order, x* are 
comoving coordinates. The photon momentum is given by = pn* with = 1. The p*'s are the 
momenta defined in the inertial frame locally defined at rest at the point (see |.22| for details). This 
choice of A in useful as it integrates out the momentum dependence from the photon one-particle 
distribution. 

In analogy with the first order calculation [29] , the full (first -|- second order) brightness equation 
can be written as: 



e/o''rf^'^(V)(l+5e(x(r,')V))i^ 

Drj 



2 ^ 

m=— 

\ m=-2 / 



Ai^^ + 4^ + IaS - ^ E V S + l^fl) y^mifi) + An . (.^(^) + 1.^(2)) 

m=-2 ^ ^ 



(19) 



where f = —neCTCL is the differential optical depth, ax is the Thomson cross section, a is the 
scale factor. A dot denotes a derivative with respect to conformal time. In the above A^*^ = 
A^'^\x{rj),h,ri), v = v(x,ri) and ^ = ^(x,rj). v is the baryon fluid velocity; Yim are the spherical 
harmonics |^ Our convention for the spherical harmonic decomposition is: 



(20) 



In writing eq. (19) we have already dropped all second order metric perturbations together with 



all products of first order terms not enhanced by 6e (the brightness equation is given in full by 



^Two technical notes. 

First, in the brightness equation in real space, we are not using the Legendre functions Pi{vb ■ n) for the 
first order terms entering the second order brightness equation as written in both [33] and [23]. Using the 
Legendre Polinomials is only justified for the first order terms in Fourier space when the choice for the z axis, 
Cz, is along the axis of symmetry, i.e. along k of the perturbation for each mode. However, at second order, 
one cannot choose a common axis of symmetry for the two Fourier modes in the convolution, and so in real 
space the expression is more complicated (though it simplifies in Fourier space) . This is why in our expression 
we are ending up with a multipole expansion in both I and m using spherical harmonics. 

Second, in solving the brightness equation, we are taking the approach of solving the full (first + second 
order) brightness equation, and then expanding the line of sight solution to second order. This approach 
gives a result where the piece of the source multiplying the (5e-containing terms in eq. Q is not tight coupling 
suppressed, and has exactly the same structure as the first order result ^ that allows us for a straightforward 
interpretation of the results. 
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eq. (76) in [22]), and the integrated Sachs- Wolf effect. Our approach is in fact to concentrate on the 
terms enhanced by 5e and to neglect metric perturbations at second order. This last point will mean 
that our calculation will be valid only for multipoles well within the horizon at recombination. 

Due to the advective derivative in the streaming part of the Boltzmann equation, the streaming 
terms are evaluated along photon trajectories, i.e. along the line of sight. Since we neglect sec- 
ond order metric perturbations and first order metric perturbations not multiplied by 6e, photon 
trajectories are straight and the position x along the line of sight is a function of r] through 

x(ri) = xq + firj . (21) 

In doing this, we neglect the effects of gravitational lensing which does not contribute significantly 
to the part of the 3-point function generated around recombination. 

Next, we can perform the line of sight integration as it is done at first order |29j using the above 
X dependence on rj. We then take the second order piece of the result, and Fourier transform it to 
find: 

rvo 

A(2)(,yo,fc,n) =2 / dr]g^°\v)S{k,7],n) . (22) 
Jo 

The visibility function is g{r]) = — fe""^, and the second order source function is given by 



S{k, r?, n) = e-^^-"(*-'') | ^A^^iv, ^) " ^ E J^^Sl^^^ mm{n) + In • ^) + 

I m=-2 



(2^)3 



5eik -k,r])- dr?'e^(^-^)-"('^'-'')f(r/')5e(r/', k - k) 
V Jvo 



X (^A^Q\k) + 4^{k) + ^A^^^\k)P2{k ■ n) + 4v^^\k) • ^) } • (23) 

The extra phase in the line-of-sight integral over Se above comes from the dependence of 5e on x{r]'), 
and not on x{r]) in ^19^ 

To obtain A'^^)(x = 0,770)) for an observer today at the origin, we integrate (22) over d^k/{2TT)^. 
We then expand the second order anisotropies in multipoles to obtain: 

^?2iVo,x = 0) = 2 j j'J dr^g'^'\r^){-z)-'^,l'^ (24) 

The source can be split into two pieces as follows: 

3 = 3" + S\ (25) 



^Ref. [27] has correctly pointed out the presence of an algebraic mistake in this point in the first version 
of this paper. After accounting for the mistake, our numerical results are marginally affected by this. 
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where 



Lm 



,(2) 



(2)f 



V2 



l,ifi 



21 + 1 



(27r)3 ' 



, (26) 



(27) 



with 



Sf{\k-k\,\k\,rj) = S'^\\k-k\,v)S^'-\\k\,r^) , 



b,2. 



(28) 



where 



S''Wk-k\,rj) = 6e{\k-k\,v)- / *d7?'f(??')<5e(|A:-A:|,7?') 



S';'\\k\,v) ^ [A^^\k)+mCk)) ^oS^fiV^^ + ^^Y'^^'hm,! - y^A^^\k)6i,2 . (30) 



(29) 



In the definition of we have already extracted the primordial density perturbation The fluctu- 



ations in 5'''* are to be meant as just the first order transfer functions. S"'^ is nothing but the first 



b.2 



order source. Also we defined fj = r] for all pieces of the source (23), except the piece containing 
the line-of-sight integral of 5e over rj' , for which f] = rf . For that piece, the exponent e"*^'^"'^^'"^*"'^) 
must be moved inside the integral over rj contained in which we do not do explicitly for clarity. 
To remind ourselves about this, we put a star after the integral sign over ry' above. 

The splitting of the source between the (a) and (6) term is quite physical. Part (a) of the source 
contains the second order monopole, dipole and quadrupole, which takes into account, as we will 
shortly see, of the perturbation to the photon diffusion scale due to bUf,. Part (5) of the source 
contains products of first order terms, and, as we will see, it corresponds to the perturbation of the 
visibility function. 



3.2 Perturbing the diffusion length: part (a) of the source 



Unlike which is the convolution of two first-order pieces, part (a) of the source (25) contains 
the lowest moments of the second order photon energy density perturbations. One should solve the 
second order Boltzmann hierarchy to obtain those exactly. In this section we will find an approximate 
solution for the first second-order multipoles, where we are able to write them as a convolution of 
two first order perturbations. The result for the source S'^ will be of the same form as 5^ in (27), 
and below we will derive the analogous expressions for S"'^ and <Sf'^. 

Here we are interested in the effect of be on the first second-order multipoles. It is easy to realize 
that this effect is tight coupling suppressed. This does not mean that the effect is negligible: even a 
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relatively small alteration of the Silk damping scale can lead, as we will see, to rather large effects. 
For this reason, the effect will be large for those high-A; modes for which the Silk damping is relevant 
^ — 0.15 Mpc""*^). Further, as we have shown in |22j . at high fe's the perturbation to 6e gets 
damped. This is due to the fact that if the mode oscillates too fast with respect to the timescales 
of recombination, the perturbation to the free electron density averages out to zero. For these two 
reasons, we can restrict our interest to the second order perturbation generated on the monopole, 
dipole and quadrupole at very high k^s as sourced by a relatively long wavelength 5e- In this limit, a 
good approximation to our solution can be found by just perturbing the diffusion damped solutions 
for the first order moments, given in e.g. [28], eq.s (8.39, 8.40). If we imagine to have a very long 
wavelength mode <5e, its main effect on the evolution of short scale perturbations will be to alter the 
local density of free electrons. We can take this into account by substituting 



ne(r?) ne(r?)(l + 5e{x, rf)) 
in eq. (8.40) of ^28j for the damping scale. We can therefore write, for / = 0, 1,2: 



A 



(i)+(2)/7: 



(fe,r/) ~ Ap^(fe,r?) il-k^ I dr]' 



1 



6(1 + i?)f (?7') 



(l + R) ' 9 



+ 



6e{x,rj') 



(31) 



(32) 



where R = 3p^^^ /{Ap^^). This is a good approximation in the limit where the wavelength of the 
temperature perturbation is much shorter than the one of the 5e perturbation. Using this equation at 
second order amounts also to neglecting the possible azimuthal dependence of A^'^\k), which again 
should be a good approximation in the same limit. Notice also that we did not need to perturb R 



in (32), since it would give rise to a second order perturbation not proportional to 5e- 



For modes with no azimuthal dependence, the relation between Legendre moments and spherical 
harmonic moments is (with k not necessarely on the z axis): 



Ai^ilr]) = y^{47r){2l + l){-ir'Ai{k,7^)Yi^M 
Using the above equation also at second order, we can write: 



A 



(l)+(2), 
Im ^ 



, rj) = A\J^{k, ■r])\l-k J dr]'kei{rj')6e{x, rj') 



where we have defined 



i\2 



+ 



{l + R{ri')) 9 



Q{l + R{i))f{i) 
In real space, this equation reads 

^™'"^'^(^.^) = AS(x,r?) + V2A;;;](x,r/) T dr?'ker(V)<^e(f , V) • 

JO 



(33) 



(34) 



(35) 



(36) 



Now, we can take into account of the mild dependence of 5e on the position by Fourier transforming 
properly the above to get 



A 



!;2+(^)(fe,r,) = A!;2(^>^) 



d^k' 
(27r)3 



dr,'\iei{r,')5e{k - k' , r,')k'^ AfJ^{k' , 77) . 



(37) 
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(2) 

We have therefore obtained an expression for which is vahd in the 'squeezed' hmit where the 
space dependence of 6e is milder than the one of a|^^ : 

AS(fc,r?) = -2 j j\^\eT{ii')5e{k-k\r^')k'\-i)-' 

X v/(47r)(2/ + l)Af)(fc', v)yL{k') (38) 

with I = 0, 1,2. 



Part (a) of the source (25) can be written as (for k along z): 



l.fh ' 

+ 25^ {vl^\k)6^,o - v^^\k)6^,i - v'^'l{k)6^,-i) + ^\2Afl\ , (39) 



where, for k = z, is defined by: 

v^'\k) = -^v^^\k)z + Y: ^'^^Ck)^ ■ (40) 
In tight coupling we have the usual relationship 

v^^Kk) = i-ir^M , (41) 

so that we can write: 



S^ik) = e-^'^<^o-v) ^(-.)y ^Yj,Jn) | U^oik, v)\,6,.,o + ^ '^U^S1(^' ^) + ^'^^,2^2^} ' 

l,rh ' 

(42) 

Though this formula is written now for generic k, we have derived it using formulas valid only for k 
along z. By using the transformation rules under rotations that we give in App. |Aj one can show 
that this expression transforms correctly under rotations, and is correct for the frame where k is 
parallel to z. Therefore it is correct in any frame. Combining this with (38) we obtain: 



= E^-'J^yir.if^)l^3Sfi\k-k\k^^^ 

l,rh ' 

where we kept fj (although it equals r] in this case), which will simplify the calculation of the 
bispecturm later on. In the above equation we defined: 

S^lk -k'\,\k'\,rj) = S'^'Wk -k'\,r,)S^^\\k'\,il) , (44) 



with 

S' 



'''\\k-k'\,T]) = - r d7]'6e{\k-k'\,T]')keT{T]') , (45) 

Jo 



S^'\\klv) ^ V4^k"A^^\k',r^)di,o + 'i^jfk'V'\k',rj)6,,, - \^^-^A^\k\^)k'H,^2 ■ 
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Thus, as anticipated, we find that even Part (a) of the source can be written in exactly the same 
form as in eq. (27). As we aheady did in the case of S(^), in we have introduced the primordial 
fluctuation The first order quantities in S""'* should be meant simply as transfer functions. Notice 
that 



S' 



a,2 



d log ko 



(46) 



where S^^^ is the first order source function. 



3.3 Line of sight solution 



(2) 

Now that we have calculated the second order source, we can put A, in a form that is keen for 



the calculation of the 3-point function. We substitute the expressions for the source, (27) and (43) 



into eq. (24) and expand the exponentials into spherical harmonics. Performing the h integral, we 
obtain: 



(%,x = o) 



(47) 



E E 



2 / T^T^ / dr^g{vV'' "'47r(2/ + 1)(2L + l)^ {21' + l){2l" + I) 



(2^)3 (27r)3 7o 



y<ju{k{m - ri))jv'{\k - fc|(% - mvnAm"n."{k - k)Y(Jk)S'^\\k - k\,v)S}'\k,v) x 

yo J y -m -m" -M J y J y M m' m J '^^ ' ' 

where i goes over a and b, and we have introduced the Wigner 3j symbols. The above result is in a 
convenient form for calculating the 3-point function. 

Before proceeding, we need to write down the first order perturbation in a similar manner. By 
analogy with the above calculation, the first order perturbation is given by: 



A2(r/o,^ = 0) = (-!)' 



(27r 



(48) 



T (1) ~* 

where A, {k,7]) is the first order transfer function, that depends only on the magnitude of k: 



A 



T,il) 



-l)V47r(2« + l) j dri ji{k{7]o - ??)) x 
k^ dr]^ 



(49) 



-re 



(1) 



-re 



4 3-point function 

We are now ready to compute the 3-point function: 

(•^iimi A/2m2 A/sma) 



\'^^hrnAhn^2^lL) + ^ Symmetric terms . 



(50) 
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Here is given by (48) and is given by (47). The main ingredients in the calculation of 
the above expectation value are the following. We assume gaussian primordial perturbations with a 
power spectrum P{k) defined as 



(51) 



Using the above expectation value in eq. (50), one obtains a term multiplied by 6^^\ki + k2 + k^)- 
The momentum dependence of 5* (with i = a, 6), is such that one can do the integral in ki trivially, 
as well as the two angular integrals on A;2,3- 

Using expression (33) for a[^, which is exact at first order, the final result of the above manip- 
ulations for the three point function arising from the source 5"* is 



l/^(2),i^(l) ^(1) ^ 

o\ hmi l2m2 hm-j,! 



E 



Ll'l 



k2dk2 f kidks 



(27r)3 J (27r) 



vo 



dr] g{r]) 



■h-V-h-l 



V47r(2/i + l) (2/2 + 1) (2/3 + 1) 



(2L + 1)(2/' + 1) V(2/ + l)(2/i + l)jv{k^{m - v))jhiHm - V)) 



xS''\k2)s}'\ks)A^^''>{k2,m)Kri''s,m)P{h)P{ks) 




X 



L I' I 
M m' rh 



-m 



r 

-m' 




(52) 



+ (2^3). 



Because of rotational symmetry, the CMB anisotropy bispectrum can be written ai 



Tjinim2inQ / \ 



where 



(7r/4)3/2 

vTOW+T 



h h h 

mi 1712 w-s 



(53) 



l2,h ' 



(54) 



is the rotationally invariant bispectrum, which is split into two pieces coming from the two pieces of 
the source function. Here we have used that the multipoles of the CMB anisotropies aim are given 
by: 



air 



1 / 47r 
4 V 2/ + 1 



d^k 
W) 



^^im{k,r]o) . 



(55) 



Using the formula for the sum over the product of two 2>j symbols given in [3F and combining (50) 



^Note that Aj^ and aim are nonlinearly related, but as usual we neglect all additional terms not enhanced 



by 5e 
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and (52), we obtain that the piece sourced by is given by 



i,mim2'ni3 



nhhh Ji+h+h 
^m\m2iny, ' 



:h-l'-l 



I' I 



h I' I 





which is a confirmation of the rotational symmetry of om' formulas. 



27r5/2 647rV(2/2 + l)(2/3 + l) Jo 

2 



(56) 



+ 5 symmetric terms 



5 Summary of the analytical results 

Let us summarize the results of our analytical treatment. It is easy to realize that our second order 



expression for A^^^ in eq. (47) exactly reproduces what we expected from eq. in the Introduction 



which implies that also the form of the bispectrum is very similar to the one given in sec. |2] In fact, 
if we start from part (6) of the source, we can see that S'^'^ is nothing but the first order source, 
while S^'^ is just equal to 6g: part (b) of the source reproduces exactly what we expected from 
the perturbation to the visibility function. Concerning the perturbation to the damping scale, this 
is given by part (a) of the source. Also in this case, S°'''^/k'^ turns out to reproduce nothing but 
the first order source (where we have not kept control of the metric perturbations, that therefore 
have disappeared). 5"'^ is just equal to ISkj^/kj^. Putting all these pieces together, and obviously 

I I (2) 

neglecting the complications coming from rotation invariance, we see that expression (47) for 
reproduces what we had anticipated in eq. 



The final expression we obtained for the bispectrum in ( 54 ) is of course quite more complicated. A 
useful simplification is obtained if we express A^ given in (49 ) using the standard photon temperature 



transfer function, Qi (given in e.g. |28j and which is returned from CMBFAST for example |29]). 
Combining 



A 



(A:,77o) = (-l)'4V47r(2Z + l)G, , 



(57) 



where T stands for transfer function, with (33), we see that a|^^ = (— i)'4Gi. Using these trans 



formations we can check that the multipoles of the CMB anisotropics defined in (55) reduce to the 
standard form (up to an irrelevant phase): 



air. 



i-i)- 



r d^k r 

J (2^7 



dnY, 



lmin)Q{k,n,r]o) 



(58) 



Further, since the source contains multipoles only up to the quadrupole, we have the following 
inequalities: < < 2 , /i — 2 < < + 2 , /a — 2 < ^3 < /a + 2 . Combining these with the 
triangle inequalities for h,l2,h, we can explicitely do the sums over l'^, I3. 
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Using (57) we find|^ 

(59) 

rvo 

grn„,.i,+i2+h / ^^^(^^ _ ^) (x)/?," (x) + (a ^ 6)] + 5 symm. , 



4 = ^/ dA;2fciei,(A:2,r?o)P(A:2)ji2(^2X^)5^'^(fc2,??o - a;) , (60) 
Pi f dhkl@i,{k3,7io)P{h) X (61) 



^J«3(^3a;)'§0^(^3,% - x) + ^ [?3ii3--l(^32;) - (^3 + l)i«3+l(^3a;)] 5'i^(^3,??0 " 2;) + 

+ ^ [3ji'3(A;3x) + i;3(A;3a;)] S'/{k3,rjo - x)^ . 

where = x for all terms in the source, except the term containing the integral over 6e in r]' , in 
which case Xfj = rjo — rj', and ji^ must be brought inside the integral over rj' in S^''^. This is our final 
expression for the bispectrum generated during recombination, which we will evaluate in the next 
section. Notice that in the limit in which the time dependence of 5e is slow compared to the width 
of the visibility function, af and a'^ can be evaluated at x = r/^ and approximately brought out of 
the X integral. What remains is proportional to the variation of Ci^ with respect to logX, where 
X are the quantities perturbed by 6e, reproducing the approximate expression for the bispectrum 
found in ( 13 ). 



6 Results and Discussion 



6.1 Signal-to- noise of the bispectrum 



We evaluate numerically eq. (59), working in synchronous gauge, using the following cosmological 
parameters: 

(Ofc, Qa, h, Tcmh,Yp, ns,T) = (0.0441, 0.742, 0.719, 2.725, 0.24, 0.963, 0.087), 

and we are now ready to compute the signal-to-noise of the bispectrum induced by perturbations to 
the recombination history. We concentrate on a full sky cosmic variance limited experiment between 
without accounting for lensing 



^min Snd /jj 

is given by (see for example |36j): 
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The Fisher matrix between two bispectra i and j 



E 

'min<'l<i2<'3<'n 



(62) 



^One may wonder why our expression for the bispectrum contains two fc-integrals, while the bispectrum 
in the local model contains 3 /c-integrals. This is due to the fact that the perturbation Se affects the CMB 
only locally in real space: this is true for the effect on the visibility function along the line of sight, and 
approximately in the squeezed limit for the effect due to the damping. Technically, this boils down to the fact 



that the second order source factorizes as in eq. ( 27 1 , which allows to eliminate one of the k integrals using 



the two delta functions coming from (51 1. 

^"^This is a good approximation up to I of order 3000 
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where A/^/j^g = 1,2,6 for triangles with no, two or three equal sides; and Ci is the CMB angular 
power spectrum. We decide to concentrate only on the bispectra from recombination and from the 
local model (with = 1). The local model is a good representative for the non-gaussianities 
one can expect from inflation (in this case multi-field inflation |3H [5^ ), and therefore it offers a 
qualitative measure of the degradation between the primordial and the cosmological signals 
The signal-to- noise ratio S/N for the recombination bispectrum is given by 

(S/Nf = Frecrec • (63) 

This assumes a negligible degradation from the local model 

We plot the total S/N{< /max) as a function of /max in Fig- 18] The minimum / we include in 
the calculation is /mm = 100. This means that, consistently with our approximation, we sum over 
modes inside the horizon at recombination, which corresponds to / ~ 70. We obtain S/N ~ 0.4 for 
Imax = 3000. This means that such a signal should not be detectable by a satellite like Planck from 
the analysis of the temperature signal alone. However, the information contained in the polarization 
signal is expected to be comparable to the one in the temperature [26], and therefore we expect 
that if our calculation were extended to include polarization, the resulting signal may be Planck 
detectable. The which in the local model generates the same signal-to-noise for the given /min 
and /max is usually referred to as and thus our result corresponds to 

/§i-«-3.5±0(l) , (64) 

where the correction 0(1) comes from what we naively expect from the terms not enhanced by 6e 
that we dropped. The effective /^l' is a convenient way to parametrize the total signal-to-noise 
but one should keep in mind that it carries no information about the shape of the bispectrum. 
The sign of is chosen according to the convention discussed below. We also compute the 
correlation coefficient between the bispectra induced by recombination and by the local model, 
rrec,ioc = Frec,iocl \J ^rel,recFf^li^^ as a function of /max- This is shown in Fig. 8 where one can see 
that for low /max it is equal to -1, and it decreases rapidly as we include higher / s. The Fisher matrix 
with the local model, for /max = 3000, is given by: 



0.19 0.011 
0.011 0.016 



(65) 



where index 1 corresponds to the non-Gaussianities from recombination, while index 2 corresponds 
to the local model. From this we can deduce that the bias of the non-Gaussianities from perturbing 
recombination into the estimate of /^l' is about 0.7 at this resolution. 

In Fig. [s] we also plot the contribution to S/N from each Imax^, and each /mm"bin by plot- 
ting d{S/Nf/dlmin and d{S/Nf/dl 

max- It is easy to see that the biggest contribution to the total 
S/N comes from squeezed triangles with 

l2-h> 2000, h ~ 200 , (66) 



Since we will find that the degradation on the local model due to recombination is relatively small, we 
can avoid to compute explicitly the degradation induced on the equilateral model, which can come from single 
field inflation. 

^^If the correlation with the local model is included in calculating the S/N , the result for /min = 100 and 
Imax = 3000 is changed only by 2%. 
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i.e. with the smah / around the first acoustic peak, and the two high Is in the diffusion damping 
tail. This result is consistent with our expectations stated at the beginning of sec. [2] based on 
considering the timescales at recombination. The fact that the signal is not dominated by triangles 
with the smallest of the Vs we consider (it decays for I < 200) is a confirmation of the viability of our 
approximation of neglecting metric perturbations. This result matches what was found explicitly in 
|21] for triangles with the smallest side out of the horizon. 




Figure 8: The plot shows the signal-to-noise for the recombination and local model bispectra; 
d{S / N)"^ / dlmax and d{S / N)"^ / dlmin (for triangles with Imin = 100 and Imax = 3000, respectively) of 
the bispectrum generated around recombination; and the correlation coefficient between the local 
model and recombination vs /max with /mm = 100. The computation was done for an ideal (cosmic 
variance limited) experiment without accounting for lensing. Clearly triangles with two high /s, and 
a low-/ around 200 dominate the recombination signal-to-noise. The local model and recombination 
bispectra start decorrelating only after high-/ modes are included. 



6.2 The shape of the non-Gaussianities 

Let us discuss the shape of the non-gaussian signal produced by fluctuations in 5e- The signal-to- 
noise per triangle, which is given by \Bi^i^i^\/ -^Ci^ Ci^ Ci^ , can be read off from eq. (63). It includes 
the additional weighting from the 3-j symbols which takes into account the phase space distribution 
of the triangles. The physical information in the bispectrum is contained in the S/N per triangle 
divided by the phase space weight of that triangle shape, which gives hi^i^i,^ / ^JGl^Cl^C\ (we prefer to 
preserve the information encoded in the sign of h). bi-^i^i^ is the reduced bispectrum and is defined as 
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^hiZ^"^^ = ^W2^'"''^''"^''^''&/ii2«3- It can be read off from QQ. The reduced bispectrum reduces 
to the bispectrum in the fiat-sky approximation [36^, and thus it contains the physical information 
of the three-point function. 

Let us comment briefly on the sign of the bispectrum. With our definition, the signs of Bi-^i^i^ and 
6i^/2^3 same. These are chosen such that they coincide with the one of the reduced bispectrum 

in [36] in the "local" model, in which a positive parameter corresponds to a negative skewness 
of the one-point distribution function of the temperature anisotropics, and to a negative reduced 
bispectrum in the Sachs- Wolfe regime (see eq. 4.24 of [36]). With our definition, a generally positive 
bispectrum (negative /nl') corresponds to high-? modes having enhanced amplitude in hot long- 
wavelength patches of the sky. 

Let us now analyze the shape of the signal. In Fig. [o] we plot /g''^ x hi-^i^i^ / {Ci^Ci^Ci^Y/"^ , for 
li ^ h (to avoid redundancy) subject to the triangle inequalities for a typical high-S/N l^. In 
the plot we choose ^3 = 1700. For comparison, we plot the analogous quantity in the case of the local 
model. The factor of l^"^ is introduced so that after normalization of the li and I2 axes by dividing 
by Z3, the integral over the square of the plotted function, multiplied by the phase space weight (see 
eq. (4.54) of [36]), is equal to the (S/N)'^ per logarithmic /max = ^3 bin. The most important feature 
to notice is that, as anticipated, the signal is peaked on squeezed triangles, but the shape is different 
from the one of the local model. 

(a) (b) 

Total Local model 




Figure 9: In panel (a) we show h^^^^ij {Ci^Ci^Ci^Y/^ for an YAgh-S/N h (we choose ^3 = 1700). 
For comparison, in panel (b) we plot the analogous quantity from the local model with 
/nl ~ T^^^ peak in the lower-right corner of panel (a) (/i ~ 200, ^2^/3 = 1700) 



corresponds to the high-S/N squeezed triangles given by (66). These triangles are such that 
they sample both the Silk damped tail, and the first acoustic peak. As for the local model, 
the signal is dominated by squeezed triangles, though the shape is clearly different. The 
acoustic oscillations can be clearly seen. The color scheme of the plots is chosen to highlight 



The bispectrum is defined up to a phase. Our definition of the bispectrum contains an extra factor of 
j^h+h+h (which is real since I1 + I2 + h is even) compared to [36j due to the extra factor of in our definition 



of aim (581. 
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For purposes of comparison with the hterature (see for example [36])) we plot in Fig. 10 the 
reduced bispectrum as a function of the smallest /, li, for isosceles configurations with the higher 
Vs, I2 = h, ranging from 2200 to ~ 2300. As /i changes, we clearly see the acoustic oscillations in 
the bispectrum. The change in value as I2 varies is also due to the acoustic oscillation. 



J2 
o 



1000 
800 
600 
400 
200 


-200 



12=13=2200,- ",2300 




10 



100 



1000 



Figure 10: The reduced bispectrum for triangles with I2 = h- The different curves correspond to 
equally spaced values of I2 = h- The bottommost (light-gray) curve is for h = h = 2200, while the 
topmost (dark-gray) curve is for h = h = 2300. The similarity between the bispectrum for these 
triangles and the temperature anisotropy power spectrum is apparent. 



In order to get some intuition on how much the signal is scale invariant, in Fig. 11 we plot the 
S/N per triangle for isosceles triangles {I2 = h)- Each curve represents a triangle with a certain 
level of squeezing, characterized by the ratio h/h, with li < I2, as we vary the overall size of the 
triangle characterized by l2- The topmost curve corresponds to the most squeezed triangles, which 
confirms that most of the signal is in that kind of triangles. However, we see that as the size of the 
triangle varies, the amplitude of the signal fluctuates a lot. This is a violation of scale invariance. 
In fact, the thick curve represents the analogous quantity for a local shape with f^f^' = 1 for the 
maximum squeezing, whose signal is generated by primordial scale invariant perturbations. We find 
that the equilateral and obtuse triangles have a completely negligible contribution. Only squeezed 
triangles with ratio /1//2 ^ 0.1 have a signal larger than the ones corresponding to Z^^' = ^• 

6.3 Physics of the bispectrum 

In sec. [2| we provided simple estimates for the origin and the size of the effects of the fluctuations 
of 6e on the bispectrum. After the full calculation, we are now able to check those results. From 
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500 1000 1500 2000 2500 3000 

Figure 11: The S/N per triangle for triangles of various shapes with I2 = h- The different curves 
correspond to different ratio of /1//2 which varies between 0.05 and 1 in equal logarithmic steps. 
Each curve terminates on the left with a dot, where h = 100, which approximately is the minimum 
I below which we can not trust our solution anymore. On the rightmost side of the plot, there are 
40 curves, and this number is reduced to the left. The bottommost (light-gray) curve corresponds 
to equilateral triangles, while the topmost (dark-gray) curve is for highly squeezed triangles with 
h/h = 0.05. The thick curve represents the analogous quantity for the local model with Z^-^' = 1 ^^'^ 
h/h = 0.05. From the plot we can see that only extremely squeezed triangles with ratio h/h ^ 0.1 
have enhanced S/N. Obtuse triangles have negligible contribution, though we do not plot them 
explicitly. 



Fig. [8] we saw that most of the signal comes from squeezed triangles with li <C ^2 — ^3- In that 
limit, the leading contributions to bi-^i^i^ come from the terms af_^{f3f^ + f3f;J and a|'^(/3;^^ -|- P^^) of 
eq. (59). These terms respectively correspond to the following combinations: low-A; and high-A; 
first order source S^^^; and low-A; 6g and high-A; S^^\ The other interesting term, high-A; 6g and low-A; 
S^^\ gives a negligible contribution as we can expect from our estimates. This is exactly the same 
structure we found in sec. [2] In particular, the fact that the highest signal from perturbations to the 
diffusion scale due to Jng, which we call 5^^, comes from squeezed triangles with a slowly varying 
Se, is a consistency check of our approximate treatment of the diffusion damping at second order in 
sec. 3.2 The contribution from the diffusion damping perturbation to the total S/N(< Imax) from 



triangles for which that approximation breaks down is ^ 10%, which is equivalent to a contribution 
of Aff^ ~ 1/2 to the bispect rum. These order one corrections to are of the order of the ones 
we expect from terms we neglect such as second order corrections from metric perturbations, and 
are therefore beyond our control. 
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Figure 12: The S/N per triangle (keeping the sign of b) for typical high-S/N squeezed triangles with 
h = h and li = 200. We also show separately the three important contributions to S/N. Notice the 
similarity between this plot and the plot obtained using the approximate treatment of the bispectrum 
in the squeezed triangle limit, Fig. [5] We reproduce all qualitative features of that plot, which serves 
as a check of our code, and allows us to pinpoint the main physical effects generating the bispectrum 
(see the text). 



In Fig. 12 we plot the three contributions discussed above to the S/N per triangle for high-5/A^ 
squeezed triangles. The bispectrum generated from 6k is generally positive and it dominates the 
one from 5g which is generally negative. Thus, the total bispectrum from recombination is generally 
positive, i.e. the amplitude of the high-/ modes is enhanced in hot long-wavelength patches of 
the sky corresponding to the first acoustic peak scale. Notice that 6g receives contributions from 
perturbations to ko due to the shift in the position of the last scattering surface {5kg in the notation 
of sec. [2]), from perturbations to the phase of the modes at the last scattering surface, and from 
perturbations to the probability for the CMB photons to originate from different times within the 
recombination era (JArea still in the notation of sec. [2]). As we anticipated, the leading contributions 
from 6kg and 5k cancel in the limit in which 6e is much slower than the timescale of recombination, 
which explains the partial cancellation between the two terms seen in Fig. 12 This can be seen also 
in Fig. [TsJ where we plot the analogous of Fig. [9] for the same three contributions. 

If we compare Fig. |5]from our approximate treatment of sec. [2] with the same one from the full 
calculation. Fig. 12, we notice that the semianalytical treatment reproduces with good accuracy the 
general trends and the oscillation frequencies, phases and relative amplitudes of the various terms. 
The additional averaging coming from the integration in time of the visibility function in the full 
calculation suppresses the bispectrum generated by (5-phase in the full calculation with respect to the 
semianalytical result, as expected. Very importantly, this tells us that we can trust both our code 
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Figure 13: The same plot as in Fig. |9]for the three interesting contributions to the recom- 
bination bispectrum. a) Bispectrum generated by low-k 5kj^; b) Bispectrum generated by 
low-A; 6kg] c) Bispectrum generated by high-/c 5Area. Each of the two main contributions is 
peaked in the squeezed hmit. Note that, as anticipated, the contribution to the bispectrum 
from the perturbations to the visibihty function is opposite in sign to the contribution from 
diffusion damping (see vertical-axis label). 



and the semianalytical analysis, and pinpoint the main physical effects generating the recombination 
bispectrum, as described in sec.|2j In particular, the low-A; 5g has a larger oscillation amplitude than 
the contribution from the low-A: 5kjy that is due to the perturbation to the phase of the first order 
source (see Fig. [6| entering in (16). The general upward trends (in magnitude) of the low-A; 
and 5g are due to perturbations to the diffusion damping generated by 5ne and 6r]r, respectively. 
By comparing the shape of the signal induced by a high-A; 5e with what we obtained in Fig. [7j we 
can confirm that the high-A; 5g bispectrum is negligible. Overall, the sum of all these contributions 
gives a bispectrum which has an ~ —3.5, which is enhanced from the first naive expectations 

of /^L ~ 1- 
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7 Summary 



The calculation of the bispectrum induced on the CMB by the standard cosmological evolution has 
become necessary to fully exploit the signal measured by the next generation CMB experiments such 
as Planck. 

In this paper we have computed in an approximate way the bispectrum due to perturbations 
in the recombination history. In a companion paper p2j we computed the first order perturbations 
to the free electron density Jg) which is one of the relevant quantities for studying the perturbed 
recombination history. This quantity does not enter in the first order calculation of the CMB 
temperature fluctuations which allows us to predict the power spectrum, but instead it enters at 
non-linear level, and therefore it can induce a bispectrum. What we found is that the perturbation 
to rif, is about a factor of 5 larger than the other perturbations, which implies the possibility for 
these fluctuations to generate a bispectrum with ~ 5, about a factor of 5 larger than naive 
expectations, and possibly detectable by next generation experiments. 

We have therefore set up to make this computation. The full second order calculation of the CMB 
bispectrum is a very hard task, with many terms that source the second order matter, radiation, 
and metric perturbations. By concentrating on modes well inside the horizon at recombination, we 
expect to be able to neglect the second order metric perturbations. This means that our results 
will be wrong by an effective ~ 0(1), which is the signal expected to be induced by metric 
perturbations Because of the large value of 6e, we have concentrated on those second order 
source terms that are proportional to 6e, and that therefore are expected to give an enhanced effect 
with respect to the rest. 

What we find is that the induced bispectrum on the temperature signal corresponds to an 
/nl — —3.5 lb 0(1). The signal is peaked on squeezed triangles with the lowest / ~ 200 and 
the highest ones at about / > 2000, and it is very far from being scale invariant. Such an 
corresponds to a signal-to-noise for an experiment like Planck of order one half. It is expected 
that the polarization signal will contain an amount of information comparable to the one in the 
temperature, and therefore such a signal may be detectable by such an experiment. 

We find that the physical origin of the signal can be understood quite simply in terms of three 
physical effects each one giving an of order a few. The first effect is due to the time shift 6r]r 
induced on the time of recombination by a low-k be mode. This induces a perturbation to the phase 
of the high-fc first order mode at the last scattering surface which is proportional to the wavenumber 
itself, and that therefore grows for more and more squeezed triangles. This growth however does not 
persist until very high A;, because there is an averaging over the width of the last scattering surface, 
and we are left with an /^l' of order a couple. 

A second effect due to the time delay is that it changes the amount of time during which photons 
can diffuse. This second effect goes together with another effect, the perturbation to the diffusion 
damping scale oc due to a low-fc be mode. In the limit in which the be perturbation is much 
slower that the timescale of recombination (which is approximately the case for low-/c modes), the 
enhanced effect from the perturbation to Ue is cancelled by the change in the time during which the 
photons diffuse, as expected. However, what is left of the perturbations is boosted in the squeezed 

Recent estimates in [37] and [3S] gives an /^l"'' ^ C(10) for equilateral configurations, which corresponds 
to roughly ^ 0(3). 
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limit by the square of ratio of the k of the high- A; mode to the damping scale ku, and again we are 
left with an effective of a few. 

By summing all of these effects, we get our final signal corresponding to ~ —3.5 it 0{1). 
Notice that this sign of corresponds to enhancing the short scale power in hotter long scale 
regions. Though physically well defined, experimentally quite large, and possibly detectable, the 
effect we find is nevertheless not much larger than the naive expectations for the effects coming 
from the terms we have neglected, which are expected to give rise to of order one. We can not 
therefore exclude that other sources might give comparable signal. In particular it is not clear to us 
that other first order quantities that are not exponentially suppressed at high k do not give rise to 
a similar effect. Therefore, in addition to providing a clear calculation of some of the second order 
effects, we think that our result further motivates the full calculation of the CMB anisotropies at 
second order. 



Note added 

While this paper was being written, a similar paper |39j appeared. Their calculation is quite similar 
to ours, though they do not include the second order monopole, dipole, and quadrupole in the 
collision term. For the bispectrum, this is equivalent to neglecting the perturbation to the 



diffusion damping due to 6ne, or in other words part (a) of our source (see our sec. 3.2). This plays 
a very important role in our calculation and in the size of the total effect. Our analytical results 
from part (6) of the source, though written in a different way, agree with theirs (see our footnote roj). 
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A Rotations 

We find it useful to summarize in this Appendix some of the transformation properties under the 
rotation group of the objects we treat in this paper. These formulas has been useful to do consistency 
checks of our calculations. 

Suppose we have a rotation R that rotates the axis by from x, y, z to x', y' , z' . A scalar transforms 

as: 

A'{k',h') = A{Rk',Rh') , (67) 

while a vector as 

if{k', n') = R-^v{Rk', Rn') . (68) 
The spherical harmonics transform according to the following relation 

Yim{Rh') = Dl^^,{R-^)Yim'{n) , (69) 
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where D^{R) is the 2Z + 1 x 2/ + 1 dimensional unitary matrix representing the rotation R. It is also 
true that 

DLoiRe.^k) = , (70) 

DL^^_^,{Rr = [-1)"' Dl^^,{R) . (71) 

Useful definitions of representations of the rotations and some symmetry relations can be found 
in do]. 
Since 

^im{k) ~ j d^nYi*^{n)A{k,h) , (72) 

we have that in a rotated frame: 

AL(^') ~ / d\'Y;^{n')A'Ck',h') ~ I d^n%:,{h')A\Rk,Rn) Dl^^,^,{R-')Air^,{Rk') . 

So we have the property under rotations: 

AUk') = D^m,m'iR-')^im'{Rk') • (73) 

With these relationships we have verified that the equations in sec. |3]and eq. (33) transform correctly. 

One can also perform a rotation of the integration variables in the Gaunt integral and find the 
following nice relationship between 3- J symbols: 

( ]=D'\ (^^±l)(*)^?'^ (i^±l)(*)z)'^ {r^')^*^ ( \ ] , (74) 

\^ mi 1712 rns J m^.m^y I m2,m2^ ' m.^,^^^ > ^ m'3 ^ ^ ' 

where means that one can even take the complex conjugate of that relation. 

One further nice usage of these rotation properties is that they allow us to compute scalar 
functions putting one of the ks on the z axis. We do not actually use this trick explicitly in the 
main part of the paper, but still, let us show this works for the case of Bi^^i^^i^{ki,k2-,k^)^ which is 
a scalar defined as 

Bh,h,hikl,k2,k^) = Y] i ] Ai^rni{ki)Ai^rn2{k2)Ai.^rn3{k3) ■ (75) 

^-^ \ mi mo m-i / 

mi,m2i'Tt3 \ / 

(Notice that we do not even need to take the expectation value to be able to make use of the 
properties under rotations, this is the reason of the tilde symbol). We have, applying a rotation 
(repeated indices are summed): 

Kh,hikuk2,ks) = ( (76) 

\ mi 1712 ^3 / 

D^^ , (R-^)D^^ , (R-^)D^' , (R~^)A' , (Rki)A' , iRk2)A' , (Rk^) . 

Then we use that 1361: 



h h ^3 \ / h h h 
nil "1-2 "13 / I m\ m'o rh'n 
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and we substitute in (76) and do the sum over mi and m2 using 



E 



nil )»7i2 



^1 h I3 

mi 1712 "23 



^1 h h 
mi m2 m3 



2/3 + 1 



obtaining: 



Bii,i2,h{ki,k2,h) 



Since by unitarity: 



7713,7713 \ ^ 

A;^^,(mi)A;^„,(flfc2)A;3„,(m3) . 



L''^ = D'h (R) , 

7713, Trig ' 7773,777,3'^ ' ' 



the sum over m3 gives a S^'^^fh'^, obtaining: 



/1 I2 h. 



rrir, rrin 



(77) 



(78) 



(79) 



Bii,i2,i3{ki,k2,h) 



E 



mi 



mo 



mo 



AU;(i?/^i)AU,(m2)A;3„,(m3) 



771, ,mo,m'. 



(80) 



i.e. the same expression in form as in (75), but computed with the rotated variables (as it should 
be for a scalar quantity). Imagine now to take the expectation value, which amounts to actually 
compute Bi^^i^^i^. In general, there will be three 3-dimensional k integral. If one chooses to put 
the second order perturbation on ki, and then chooses the reference frame where Rki = kiCz, then 
one can use the analytical results found in the frame where k is along z. Then, one integrates over 
Rk2,Rk3, and finally the integration on ki is trivial, giving just An. With this technique, we have 



been able to reproduce the result of eq. ( 56 ) 
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